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ABSTRACT 



We describe the goals, design, implementation, and initial progress 
^ ■ of the UKIRT Infrared Deep Sky Survey (UKIDSS), a seven year sky 

Qh' survey which began in May 2005. UKIDSS is being carried out using the 

UKIRT Wide Field Camera (WFCAM), which has the largest etendue of 
\^ • any IR astronomical instrument to date. It is a portfolio of five survey 

^ \ components covering various combinations of the filter set ZYJHK and 

C^ . H2. The Large Area Survey, the Galactic Clusters Survey, and the Galac- 

tic Plane Survey cover approximately 7000 square degrees to a depth of 
K~18; the Deep Extragalactic Survey covers 35 square degrees to K'^21, 
and the Ultra Deep Survey covers 0.77 square degrees to K~23. Summed 
^ ' together UKIDSS is 12 times larger in effective volume than the 2MASS 

survey. The prime aim of UKIDSS is to provide a long term astronom- 
ical legacy database; the design is however driven by a series of specific 
goals - for example to find the nearest and faintest sub-stellar objects; to 
discover Population II brown dwarfs, if they exist; to determine the sub- 
stellar mass function; to break the z=7 quasar barrier; to determine the 
epoch of re-ionisation; to measure the growth of structure from z=3 to the 
present day; to determine the epoch of spheroid formation; and to map 
the Milky Way through the dust, to several kpc. The survey data are be- 
ing uniformly processed. Images and catalogues are being made available 
through a fully queryable user interface - the WFCAM Science Archive 
(WSA : http://surveys.roe.ac.uk/wsa). The data are being released in 
stages. The data arc immediately public to astronomers in all ESO mem- 
ber states, and available to the world after eighteen months. Before the 
formal survey began, UKIRT and the UKIDSS consortium collaborated in 
obtaining and analysing a series of small science verification (SV) projects 
o 9nn7 R A c; ^° complete the commissioning of the camera. We show some results from 

these SV projects in order to demonstrate the likely power of the eventual 
complete survey. Finally, using the data from the First Data Release we 
assess how well UKIDSS is meeting its design targets so far. 

Key words: surveys, infrared: general 
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1 INTRODUCTION 

The UKIRT Infrared Deep Sky Survey (UKIDSS) can 
be considered the near-infrared counterpart of the 
Sloan Digital Sky Survey (SDSS; York et al. 2000). 
It does not cover the whole sky, but is many times 
deeper than the Two Micron All Sky Survey (2MASS; 
Skrutskie et al. 2006). It is in fact not a single sur- 
vey but a survey programme combining a set of five 
survey components of complementary combinations 
of depth and area, covering several thousand square 
degrees to A' ~ 18, 35 square degrees to A" ~ 21, and 
0.77 square degrees to if ~ 23. The survey uses the 
Wide Field Camera (WFCAM) on the 3.8m United 
Kingdom Infra-red Telescope (UKIRT). WFCAM has 
an instantaneous field of view of 0.21 square degrees, 
considerably larger than any previous IR camera on 
a 4m class telescope, along with a pixel size of 0.4". 
The tip-tilt system on UKIRT delivers close to natu- 
ral seeing (median size 0.6") across the whole field of 
view. The combination of large telescope, large field 
of view, and good image quality, makes possible a sur- 
vey of considerably greater scope than 2MASS. The 
various surveys employ up to five filters ZYJHK cov- 
ering the wavelength range 0.83 — 2.37 fim and extend 
over both high and low Galactic latitudes. The survey 
began on 2005 May 13, and is expected to take seven 
years to complete. The survey is being carried out by 
a private consortium but is fully public with no pro- 
prietary rights for the consortium. Data products are 
being being released in stages, with the intention of 
having these roughly twice a year. 



1.1 Origins and nature of project 

The UKIDSS survey concept first emerged in 1998 
while making the funding case for the WFCAM in- 
strument itself, but eventually became a formal ref- 
ereed proposal to the UKIRT Board in March 2001, 
submitted by a consortium of 61 UK astronomers. 
This included a commitment to making data avail- 
able immediately to all UK astronomers (not just con- 
sortium members), plus specified individual Japanese 
consortium members, and available to the world after 
a year or two. (See section [V.4I for the final data re- 
lease policy.) Later, during the UK's entry in to ESQ, 
it was agreed that astronomers in all ESQ member 
states would have the same data rights as UK as- 
tronomers, and at the same time, membership of the 
consortium was extended to any interested European 
astronomers. Consortium membership now stands at 
130. Note that individual astronomers are members, 
not their institutions. 

The project is unusual compared to previous 
large survey projects, being neither private, nor con- 
ducted by a public body on behalf of the commu- 
nity. UKIDSS relies on the separate existence of three 
things, (i) The UKIRT observatory, operated as part 
of the UK's Joint Astronomy Centre (JAC). (n) The 



WFCAM instrument, built at the Astronomy Tech- 
nology Centre (ATC) at the Royal Observatory Ed- 
inburgh, as a funded PPARC project. Note that WF- 
CAM was built as a common user instrument to be 
part of the UKIRT suite of instruments. The UKIDSS 
consortium is essentially the largest single user, (iii) 
The pipeline and archive development project, run by 
the Cambridge Astronomy Survey Unit (CASU) and 
the Edinburgh Wide Field Astronomy Unit (WFAU), 
and funded by several different PPARC grants. This 
data processing development is part of the VISTA 
Data Flow System (VDFS) project, with the WF- 
CAM pipeline and archive being seen as an interme- 
diate step. Note this data processing project deals 
with all WFCAM data, not just the UKIDSS data. 

The consortium has no data privileges, and does 
not plan the scientific exploitation of the survey; its 
purpose is to make the survey happen, on behalf of 
the ESQ communuity. The aim of the UKIDSS con- 
sortium is then to produce the scientific design for 
the survey; to win the telescope time necessary; to 
plan the implementation of the survey, liaising with 
the other bodies above; to staff the observing imple- 
mentation; to define the necessary Quality Control 
(QC) filtering stages to produce final survey prod- 
ucts; to assist the data processing team as necessary 
in producing stacked and merged survey products; 
and finally to document the production of the survey 
data in scientific publications and other technical pa- 
pers. A number of individuals in ATC, JAC, CASU 
and WFAU are also members of UKIDSS, so that 
the liaison with the camera construction and data 
processing projects, as well as telescope operations, 
has been well motivated. A clear relationship with 
UKIDSS has been built into each of these projects. 
For example the science requirements document for 
the pipeline and archive emerged from consultation 
with UKIDSS; and the commissioning schedule for 
WFCAM include a "science verification (SV)" phase 
following standard tests, in which the survey imple- 
mentation described in section |3] could be tested and 
refined. 



1.2 Technical reference papers 

This paper is one of a set of five which provide the ref- 
erence technical documentation for UKIDSS. It sum- 
marises the scope, goals, and overall design of the 
survey, along with a brief discussion of implementa- 
tion methods, progress to date, and presentation of 
early "science verification" data. The other four pa- 
pers, described briefly below, are Casali et al. (2007), 
Hewett et al. (2006), Irwin et al. (in preparation) and 
Hambly et al. (in preparation). In addition to these 
five core reference papers, each data release will be 
accompanied by a paper detailing its contents and im- 
plementation information. The first two of these are 
Dye et al. (2006) for the "Early Data Release (EDR)" , 
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and Warren et al. (2007) for the "First Data Release" 
(DRl). 

Casali et al. (2007) describe the survey instru- 
ment, WFCAM. A short summary is given in sec- 
tion |3l At the time of commissioning, 2004 Novem- 
ber, the instrument etendu^ of 2.38 m^ deg^ was the 
largest of any near-infrared imager in the world. 
The Canada France Hawaii Telescope WIRCam in- 
strument (Puget et al. 2004) covers a solid an- 
gle of 0.1 deg^ per exposure giving an etendue of 
l.llm^deg^. WFCAM is likely to remain as the 
near-infrared imager with the largest etendue in the 
world until completion of the near-infrared camera 
for VISTA (Dahon et al. 2004). 

The data flow system for WFCAM is described 
by Irwin et al. (in preparation) and Hambly et al. (in 
preparation) . A summary is given in section [S] The 
very high data rate (1TB per week) requires a highly 
automated processing system that removes instru- 
mental signature, produces object catalogues, and in- 
gests into a fully queryable WFCAM Science Archive 
(WSA). It is expected that nearly all science analysis 
of UKIDSS wiU be initiated through the WSA. 

The photometric system is described in Hewett 
et al. (2006). The survey uses five broadband filters, 
ZYJHK. The JHK passbands are as close as possi- 
ble to the MKO system; the Z passband is similar to 
the SDSS z' passband, but has a cleaner red tail. The 
Y passband is a new one centred at 0.97 /im, which 
bridges the gap between Z and J. The exact passband 
was designed with the aim of discriminating between 
high-redshift quasars and brown dwarfs. 

Hewett et al. present the measured passband 
transmisssions, and use synthetic colours of various 
classes of astronomical object to produce expected 
colour equations between certain WFCAM, SDSS, 
and 2MASS filters. Table [T] provides summary in- 
formation on the ZYJHK passbands. A later paper 
(Hodgkin et al. in preparation) will report on the pho- 
tometric calibration of the UKIDSS survey and colour 
equations determined on the sky from standard star 
observations. 



1.3 Plan of paper 

This paper begins with a description of the science 
goals of UKIDSS, and some illustrations of how the 
survey design will achieve them. We then describe the 
practical implementation of the survey - the tiling and 
jittering patterns, exposure times, calibration plan, 
and so on, in the context of the camera properties 
and the UKIRT operating procedures and software. 
Next we consider the detailed design of the individual 
survey components - areas, field selection, filters and 
scheduling. We also describe the staging of UKIDSS 



^ product of telescope collecting area, and solid angle of 
instrument field of view, sometimes called grasp 



Table 1. WFCAM passbands used for UKIDSS. Columns 
show the name, the effective wavelength, (Schneider, 
Gunn & Hoessel, 1983), the short and long wavelength 
limits, and the width, i.e. the difference between these 
two wavelengths. The wavelength limits are set by the to- 
tal system transmission, and defined as 50% transmission 
relative to the peak. Complete transmission functions are 
shown in Hewett et al (2006). 



band 


Kti 


range 


width 




^m 


/^m 


^m 


Z 


0.8817 


0.836-0.929 


0.093 


Y 


1.0305 


0.979-1.081 


0.102 


J 


1.2483 


1.169-1.328 


0.159 


H 


1.6313 


1.492-1.784 


0.292 


K 


2.2010 


2.029-2.380 


0.351 



in a two year plan and final seven year plan. We then 
summarise the data processing arrangements, which 
as described above are pursued as a formally separate 
project, but which of course are crucial to the scien- 
tific success of UKIDSS. Following this we present 
some example data and simple analysis from the sci- 
ence verification phase of UKIDSS, and point towards 
the expected final data quality. Finally we summarise 
progress to date, including an analysis of how well 
UKIDSS is meeting its design targets, describe the 
plan for publication of the data, and provide links to 
more detailed information about UKIDSS, WFCAM, 
and the science archive. 



2 SCIENCE GOALS 

2.1 General goals 

The primary goal of UKIDSS is to produce IR sky 
atlases and catalogs as a fundamental resource of 
lasting significance analogous to the various Schmidt 
photographic sky surveys of the 1970s and onwards 
(Hambly et al. 2001 and references therein), and to 
the SDSS survey of modern times (York et al. 2000). 
None of our survey components covers the whole sky, 
but none the less each component deserves the term 
"atlas", as the volume surveyed, and the number of 
objects detected are comparable to the above optical 
surveys, and each survey maps out some significant 
part of the universe - the solar neighbourhood, the 
Milky Way, the local extragalactic universe, the uni- 
verse at z=l, and the universe at z=3. Each of the 
component surveys is many times larger than any ex- 
isting IR survey at comparable depth. 

The strength of a survey is of course its poten- 
tial for multiple use over many years, but this general 
aim does not fix the best combination of area, depth, 
and wavelength coverage. In a Euclidean volume, for 
a given total time, a shallow survey always produces 
a larger sample size than a deep one; but specific sci- 
ence goals often require a given depth, for example to 
detect galaxies at a given redshift; and for relatively 
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Survey Name 

Notes 



Large Area Survey (LAS) 
Matched to SDSS areas 
Deep Extragalactic Survey (DXS) 
Four multiwavelength fields 
Ultra Deep Survey (UDS) 
Subaru/XMM Deep Survey Field 
Galactic Plane Survey (GPS) 
Northern Plane, \b\ < 5° 
Galactic Clusters Survey (GCS) 
Ten large open clusters 



Passbands Depth Area 

K mag sq.deg 



YJHK 

JK 

JHK 

JHK 

ZYJHK 



18.2 



20.8 



22.8 



18.8 



18.6 



4028 



35 



0.77 



1868 



1067 



Seeing Phot Astrom Ellipticity 

FWHM RMS RMS 

arcscc mag arcsec 



<1.2 


<0.02 


<0.1 


<0.25 


<1.3 


<0.02 


<0.1 


<0.25 


<0.8 


<0.02 


<0.1 


<0.25 


<1.0 


<0.02 


<0.1 


<0.25 



<1.2 



<0.02 



<0.1 



<0.25 



Table 2. General design of each component of UKIDSS. These are the planned survey goals as originally approved by the 
UKIRT Board, together with survey quality parameters applied during preparation for the First Data Release. The final 
areas and depths will depend on a variety of factors including survey efficiency and UKIRT time allocation. 



deep surveys, neither the Milky Way nor the universe 
at large are Euclidean volumes. As we cannot predict 
all future uses of the UKIDSS survey databases, the 
general idea is to pursue a "wedding cake" strategy, 
dividing the time between a large shallow survey, a 
medium sized fairly deep survey, and a small very 
deep survey, and including targeted observations of 
the Galactic Plane and nearby clusters. 

The general features of the five planned survey 
components are described in Table [2] (More details 
are given in section |4] and in the data release papers. 
Dye et al. (2006) and Warren et al. (2007a,b)). The 
intention is to achieve uniformity of depth across each 
survey. |j Separate limits apply in each band of course. 
The photometric and astrometric accuracy goals are 
requirements set by the UKIDSS consortium on the 
VDFS data processing system. They are maximum 
absolute errors, and so include any non-uniformity 
across the survey. Seeing and ellipticity 'goals' in the 
table are a simplified version of the cuts that have 
been applied in quality control filtering during prepa- 
ration for the first data release; the actual cuts vary 
between bands, and the great majority of dataframes 
easily pass these cuts. 

In the following subsections we summarise the 
scientific goals of each survey component, and overall 
survey quality goals. In section [^ we describe the de- 
sign of each survey component - areas, depths, field 
locations, filters, implementation strategy - needed 
to achieve those goals. In section [7] we summarise the 
progress so far towards the design goals, and the ac- 
tual achieved quality. 

The scope of the surveys is illustrated in an 
interesting way in Fig. [l] For sky limited observa- 
tions in the K-band, one can show that for a uni- 
form source population in Euclidean space, the ef- 
fective volume surveyed is proportional to the quan- 

^ The depths quoted in this table and elsewhere are mag- 
nitudes in the Vega system for a point source which attains 
a 5(T detection within a 2"aperture. 



tity Area x IQ"''" . On the other hand, the time to 
reach depth K is proportional to lO"'*^. So, as one 
surveys the same area deeper, the effective volume 
increases slowly, but the overall information gained 
(for example on colours of brighter sources) increases 
somewhat faster. For fixed observing conditions, the 
quantity 10°*^ is proportional to etendue x instru- 
ment throughput X time spent, which we can think 
of as overall 'effective survey resource'. The largest 
existing multiband near-IR survey in terms of both 
quantities is 2MASS. In Fig. [l] we have normalised 
the computed values for each of the five UKIDSS ele- 
ments, to the 2MASS values using the 2MASS limit of 
K=15.50 (Skrutskie et al. 2006). Viewed in this way, 
each of the 5 K-band surveys is between 10 and 30 
times larger than 2MASS in terms of resource, and, 
except for the UDS, is a few times larger in volume. 
Summed over the whole programme, UKIDSS is 70 
times larger than 2MASS in terms of resource, and 
12 times larger in terms of volume. 

In practice, achieving the science goals of 
UKIDSS will depend not just on the depths and areas 
achieved, but also on the survey quality - complete- 
ness, reliability, number of artefacts and so on. These 
issues are discussed in section [S] with some prelimi- 
nary results given in section [7] 



2.2 Headline science goals 

To develop a plan beyond this very general concept, 
we encouraged the formation of groupings within the 
consortium to promote distinct survey components 
and science goals. The whole consortium then de- 
bated the science goals and designs proposed by these 
groups, looked for overlaps, and made compromises 
until we felt we had a balanced strategy. (The re- 
sulting survey component designs are presented in 
Section ID) As part of this scientific debate, we also 
agreed headline science goals, which were used to 
drive the specific designs of the survey components. 
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Figure 1. Illustration of the scope of the five UKIDSS 
survey components, and their sum, by comparison with 
2MASS. The comparison is made in terms of effective 
volume, and effective resource, for the K band, computed 
as described in the text 



The most important of these specific goals are as fol- 
lows : 

• to find the nearest and faintest sub-stellar ob- 
jects 

• to discover Pop II brown dwarfs, if they exist 

• to determine the substellar mass function 

• to break the z=7 quasar barrier 

• to determine the epoch of re-ionisation 

• to construct a galaxy catalogue at z=l as large 
as the SDSS catalogue 

• to measure the growth of structure and bias from 
z=3 to the present day 

• to determine the epoch of spheroid formation 

• to clarify the relationship between quasars, 
ULIRGs, and galaxy formation 

• to map the Milky Way through the dust, to sev- 
eral kpc 

• to increase the number of known Young Stel- 
lar Objects by an order of magnitude, including rare 
types such as FU Orionis stars 



in the J-band, with an average epoch separation of 
3.5 years, and a minimum separation of 2 years, to 
detect proper motions of low mass objects and thus 
their kinematic distances. The final J depth target is 
therefore J=19.9. (The final achieved depth may of 
course be different; this is discussed in [7]) 

The Large Area Survey, when combined with the 
matching SDSS data, will produce a catalogue of a 
half a million galaxies with colours and spectra, and 
several million galaxies with photometric redshifts; 
will detect thousands of rich clusters out to z=l; will 
find ten times more brown dwarfs than 2MASS, will 
probe to much fainter objects, and can get statistical 
ages and masses from kinematics; and will produce a 
complete sample of 10,000 bright quasars, including 
reddened quasars, using the K excess method (War- 
ren, Hewett and Foltz 2000). 

We are particularly driven however by three es- 
pecially exciting prospects, (i) A search for the the 
nearest and smallest objects in the solar neighbour- 
hood. The LAS is deep enough to detect brown dwarfs 
and young (< 5 Myr) free floating planets with as lit- 
tle as 5 Jupiter masses out to distances of tens of par- 
sees. The LAS should find brown dwarfs even cooler 
than T dwarfs, Te <700K, a new spectral class tenta- 
tively named Y dwarfs (Leggett et al. 2005). (ii) The 
combination of IR and optical colours, and large ex- 
pected proper motions, will allow the LAS to find halo 
brown dwarfs if they exist, testing the universality of 
star formation processes, and the formation history 
of the Milky Way. (Calculations based on Burrows et 
al. (2003) indicate that halo brown dwarfs will have 
an absolute Y magnitude >16 and so could be found 
to distances of 60pc.) (iii) We hope to find quasars at 
z = 7 and to detect the epoch of re-ionisation. SDSS 
have found z=5-6 quasars by "i' drop-out". Beyond 
z=6 quasars become rapidly redder, indistinguishable 
from brown dwarfs in standard colours, and too faint 
to be in the SDSS z' survey. We therefore intend to 
undertake a survey in the new Y filter to match our 
JHK survey. Extrapolating the evolution function of 
Fan et al. 2001 to predict quasar numbers detected in 
both Y and J to 10 a, the LAS should find 10 quasars 
in the range z=6-7 and 4 in the range z=7-8. Figure[2] 
illustrates how the UKIDSS filter set can distinguish 
cleanly between very cool brown dwarfs and very high 
redshift quasars. 



2.3 Goals of the Large Area Survey (LAS) 

The Large Area Survey (LAS) aims to map as large a 
fraction of the Northern Sky as feasible (4000 square 
degrees) within a few hundred nights, which when 
combined with the SDSS, produces an atlas cover- 
ing almost an order of magnitude in wavelength. Fur- 
thermore a huge number of objects will already have 
spectroscopic data from the SDSS project. The tar- 
get depths of the basic shallow survey are Y=20.3, 
J=19.5, 11=18.6, K=18.2. We also plan a second pass 



2.4 Goals of the Galactic Plane Survey 
(GPS) 

The Galactic Plane Survey aims to map half of the 
Milky Way to within a latitude of ±5° . Given the dec- 
lination constraints of UKIRT, we can survey 1=15°- 
107° and l=141°-230°. Owing to interest in recent re- 
sults from multi- waveband observations of the Galac- 
tic Centre region (eg. Wang et al. 2002; Hasegawa et 
al. 1998) the survey region has been extended south 
to include the l=-2° to 15° region in a narrow strip 



© 2007 RAS, MNRAS 379, 1599-1617 



6 Lawrence et al. 



:< 1.5 


1 1 1 


1 1 1 1 


1 1 1 1 1 1 1 


^ 


. 




: 


X 


- 




I f\ J 


2 1 


~ /V^ 






j^ 


r^ 






V 


'■- 




/ V_ rT, /*""2^ V 


.£ 0.5 


- 


' \/ , \ 


>'' 'i l~7 ^^ — ih 


relat 


:,//, 


t^k.w/..^ 


6 


0.8 1 1.2 






wavelength (/xm) 



Figure 2. P/ot illustrating the usefulness of the Y band 
for finding cool brown dwarfs and quasars of very-high 
redshift (z > 6.4J. Filter curves are total system through- 
put (above atmosphere to detector), normalised to the 
peak, for the SDSS i' and z' bands (from Fan et al. 2001), 
and the WFCAM Y and J bands (from Hewett et al. 
2006). The dashed curve is the spectrum of the T6 brown 
dwarf SDSS J162414. 37-1-002915. 6 (from Leggett et al. 
2000), and the solid curve is a model spectrum of a quasar 
at z=7. High-redshift quasars and brown dwarfs may be 
identified by the very sharp spectral discontinuity in mov- 
ing from the optical (i', z') to the near-infrared (Y, J), 
while the quasars may be distinguished from the brown 
dwarfs, because they are somewhat bluer in Y-J colour. 



at 6 = ±2°. The target depths are J=19.9, H=19.0, 
K=18.8, with the K depth made up of three separate 
epochs with a depth of K=18.2 each time|j This is 
deep enough to probe the IMF down to M~ 0.05 M0 
in star formation regions within 2 kpc of the sun, to 
detect stars below the main sequence turn off in the 
galactic bulge, and to detect luminous objects such as 
OB stars and post-AGB stars across the whole galaxy. 
The K band repeats are for both extra depth and to 
detect highly variable objects and locate nearby ob- 
jects through their proper motions. In addition we 
will make a three epoch narrow band II2 survey in a 
300 square degree area of the Taurus- Auriga-Perseus 
molecular cloud complex (with JHK data also). This 
survey area closely follows the region of molecular 
emission detected by Ungerechts & Thaddeus (1987). 
Like the high latitude LAS, the GPS has its 
prime importance as a fundamental resource for fu- 
ture astronomy. The survey depths are close to being 
confusion limited, so this survey is unlikely to be su- 
perseded until a high resolution wide angle camera is 
placed in space. We expect to detect several times 10^ 
sources in total. However, there are a number of im- 
mediately expected results, which will be achieved in 
combination with data from multiwaveband galactic 
surveys from many facilities. There will be particular 
benefit from surveys planned or in progress with the 



•^ These depths refer to uncrowded regions well away from 
the Galactic Centre and a few degrees out of the plane. In 
crowded regions the survey will be less deep, due to added 
background noise from unresolved stars (see section 6.2) 



Isaac Newton Telescope (optical), SPITZER space 
telescope (infrared), CHANDRA and XMM- Newton 
(X-ray), the VLA (radio, especially the 5 GHz COR- 
NISH survey), HERSCHEL and SCUBA-2 (submm) 
and AKARI (far infared). The following list illus- 
trates some of the expected results. (1) An increase in 
the number of known Young Stellar Objects (YSOs) 
by an order of magnitude and measurement of the 
duration of the YSO phase as a function of mass 
and environment. (2) Star formation regions will be 
mapped throughout the Milky Way, measuring the 
environmental dependence of the IMF to low masses 
and estimating the overall star formation rate of the 
galaxy. (3) Rare or brief duration variables will be 
found in significant numbers, aiding the study of phe- 
nomena such as FU Orionis variables. Luminous Blue 
Variables and unstable post-AGB stars undergoing 
thermonuclear pulsations. (4) Thousands of evolved 
objects such as protoplanetary nebulae and planetary 
nebulae will be found, a huge increase over previous 
samples. This will be achieved by selecting candidates 
from far IR surveys such as the Akari (ASTRO-F) 
survey. New detections are expected to be commonest 
for small protoplanetary nebulae and very large low 
surface brightness planetary nebulae, which would re- 
spectively be unresolved or undetected by 2MASS. 
(5) Many stellar populations will be mapped to large 
distances through the Milky Way extinction, measur- 
ing the scale height versus stellar type and mapping 
poorly measured regions of the arms and warp. (6) 
The IR counterparts of hundreds of X-ray binaries, 
thousands of CVs, and thousands of coronally active 
stars will be identified and source lists provided for 
regions yet to be mapped by X-ray satellites. 

2.5 Goals of the Galactic Clusters Survey 
(GCS) 

The Galactic Clusters Survey (GCS) aims to survey 
ten large open star clusters and star formation as- 
sociations, covering a total of 1067 sq.deg. using the 
standard single pass depth (see section IT2|l plus a sec- 
ond pass in K for proper motions, giving a depth of 
Z=20.4, Y=20.3, J=19.5, H=18.6, K=18.6. The tar- 
gets are all relatively nearby, are at intermediate to 
low Galactic latitudes and are several degrees across. 
The GCS is the most targeted of our surveys, 
being aimed at the crucial question of the sub-stellar 
initial mass function (IMF). Our current knowledge 
of the IMF is illustrated in Fig. [S] by comparing 
measurements of the Pleiades with various functional 
forms. (Note that determining whether the Pleiades 
is typical or not is part of the goals of the survey). 
The stellar IMF is well determined down to the brown 
dwarf boundary but is much less well known below, 
and it is not known whether the IMF as a whole is 
universal or not (the current state of research into ul- 
tra low-mass star formation is described in Martin & 
Magazzii 2007). The mass limit reached varies some- 
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Figure 3. Various extrapolations (dotted lines) of the 
Pleiades mass function (after Hambly et al. 1999) illus- 
trating uncertainties in the behaviour of the Mass Func- 
tion (MF) in the brown dwarf (BD) regime. (Vertical axis 
is cluster stars per unit mass.) The most recent surveys 
have probed the mass range 0.01 to 0.1 solar masses, us- 
ing a variety of techniques, and have produced a range 
of different forms of the MF using heterogeneous datasets 
with varying degrees of completeness. The GCS aims to 
obtain maximal completeness in ten targets to settle the 
questions as to the form and universality (or otherwise) 
of the MF in the BD regime. 



what from cluster to cluster, but is typically around 
Ml ~ 30Mj, where Mj is the mass of Jupiter. The 
number of objects expected to be detected in the 
range Ml - Ml + IQMj ranges from 100 to 3000 
for the range of possible mass function models, show- 
ing how well we will constrain the IMF compared to 
current knowledge. 

To find extreme objects — the very nearest ex- 
amples, the lowest mass objects — the Large Area 
Survey is better. But to measure the substellar IMF, 
one wants to target the 30 - 100 Mj region, and to 
obtain masses one needs both a distance and an age, 
for which mapping clusters is ideal. This approach 
has of course already been started (e.g. Moraux et a7. 
2007 and references therein). Our survey improves on 
current studies not by going deeper but by collecting 
much larger numbers, and examining objects formed 
in environments having a range of ages and metallic- 
ities, to examine the question of universality. 



2.6 Goals of the Deep Extragalactic Survey 
(DXS) 

The Deep Extragalactic Survey (DXS) aims to map 
35 sq.deg. of sky to a 5-cr point-source sensitivity of 
J=22.3 and K=20.8 in four carefully selected, multi- 
wavelength survey areas. Central regions of each field 
will also be mapped to H=21.8. The primary aim of 
the survey is to produce a photometric galaxy sample 
at a redshift of 1-2, within a volume comparable to 



that of the SDSS, selected in the same passband (rest 
frame optical). Figure|4]shows measured K magnitude 
versus redshift for galaxies in the Hawaii Deep Fields 
(L.Cowie, private communication). This shows that 
to achieve a sample such that the median redshift is 
z ~ 1 requires measuring galaxies with i^ ~ 20 and 
so going to a point source depth of K~21. Such a 
sample will allow a direct test of the evolution of the 
galaxy population and determine how galaxies of dif- 
ferent types (passive, star-forming, AGN) trace large 
scale structure (their bias). Each of these properties 
can be predicted from cosmological simulations so the 
DXS will set tight constraints on these models in vol- 
umes less susceptible to cosmic variance than previ- 
ous, narrow-angle surveys at this redshift. The sam- 
ple will also enable the selection of clusters of galaxies 
in this redshift range, where cosmological models pre- 
dict numbers to be sensitive to the total mass density 
of the Universe, f2o- 

The number of deep, multi- wavelength survey 
fields has increased dramatically in the past 5 years 
with the up-grade of existing facilites (e.g. Megacam 
on CFHT and VIMOS on the VLT) and new satellite 
missions (e.g. Spitzer, GALEX, XMM-Newton and 
Chandra). Each of these facilities has current surveys 
of 2-40 sq.deg. of contiguous area to levels where 
many of the counterparts are intrinsically faint in 
the optical (R>24) due to a combination of redshift 
and/or intrinsic dust obscuration, but are relatively 
red (R-K>4). Therefore deep NIR imaging (K~21) 
over tens of square degrees is required to fully char- 
acterise these dusty and/or distant objects. Looking 
ahead to the end of this decade and the completion 
of UKIDSS, there will be many more complementary 
surveys on these scales in other wavebands such as 
the far- infrared (Herschel), sub-mm (SCUBA-2), ra- 
dio (EVLA and eMerlin) and Sunyaev-Zel'dovich (SZ 
telescope and AMI). The legacy potential of the DXS 
was a key driver for the science case and the field 
selection (see section 4.4). 

2.7 Goals of the Ultra Deep Survey (UDS) 

The Ultra Deep Survey (UDS) aims to map 0.77 sq. 
degrees to a Scr point-source sensitivity of J=24.8, 
11=23.8, K=22.8. Such depths are required to reach 
typical L* galaxies at z = 3. Covering an area one 
hundred times larger than any previous survey to 
these depths, this will provide the first large-volume 
map of the high-redshift Universe (30 x 30Mpc by 2 
Gpc deep at 2 < 2 < 4). 

Deep near-infrared surveys are crucial for obtain- 
ing a more complete census of the Universe at these 
epochs. In particular, galaxies which are reddened by 
dust or those which appear red due to old stellar pop- 
ulations may be completely missed by standard opti- 
cal surveys. From the UDS we anticipate over 10,000 
galaxies at z > 2, allowing detailed studies of the lu- 
minosity functions, clustering and multi-wavelength 
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16 17 18 

K MAGNITUDE 

Figure 4. The redshift distribution of a K selected galaxy 
sample from the Hawaii Deep Fields. This is an updated 
version of the figure in Songaila et al (1994), kindly pro- 
vided by L.Cowie and collaborators. The solid symbols 
show spectroscopic redshifts from the Hawaii Deep Fields, 
which have been completely observed to K=20 (dashed 
line) though only identified objects are shown. The open 
diamonds show the spectroscopically identified objects in 
the Hubble Deep Field, while the crosses show all the re- 
maining objects at their photometric redshift. The solid 
line shows the K magnitude of a 2L* unevolving Sb galaxy. 



SEDs over a large, representative volume. A major 
goal, together with the DXS and local surveys, is to 
measure clustering as a function of stellar mass and 
cosmic time, which will provide very powerful tests of 
models for biased galaxy formation and the growth of 
structure. 

The UDS is also designed to address one of 
the major unsolved problems in modern astronomy, 
which is to understand when the massive elliptical 
galaxies are formed. A key test will be to determine 
the co-moving number density of the most massive 
galaxies at various epochs, particularly at z > 2. 
This requires a combination of both depth and area 
which has previously been impossible to achieve. If 
the density of massive galaxies (more massive than 
local L* ellipticals) is similar to that of today, we 
should see ~ 1000 per square degree. Current semi- 
analytic models predict an order of magnitude fewer. 
Our goal is to directly measure the build-up of this 
population over cosmic time. 

The survey field chosen for the UDS is the Sub- 
aru/XMM Deep Field, which has a wide range of 
multiwavelengh data available, including deep radio 
observations from the VLA, submm mapping from 
SCUBA, mid-IR photometry from Spitzer, deep opti- 
cal imaging from Subaru Suprimecam and deep X-ray 
observations from XMM-Newton. When combined, 
these will enable detailed studies of the relationship 
between black hole activity, dust-dominated ULIRGs 
and IR-selected massive galaxies using an unprece- 
dented high-redshift sample. 



3 IMPLEMENTATION WITH THE 
UKIRT WIDE FIELD CAMERA 

3.1 General characteristics of telescope and 
camera 

UKIDSS is implemented using the Wide Field Cam- 
era (WFCAM) on the United Kingdom Infrared Tele- 
scope (UKIRT), which is operated by the Joint As- 
tronomy Centre (JAC), an establishment of the UK's 
Particle Physics and Astronomy Research Council 
(PPARC). General technical details for UKIRT are 
given on the JAC websitqj. It is an infra-red ded- 
icated 3.8m telescope operating at the summit of 
Mauna Kea in Hawaii. Of particular importance is a 
tip-tilt secondary, which primarily removes dome and 
windshake effects on seeing, delivering close to free- 
atmosphere seeing (half arcsecond on many nights) 
across the whole WFCAM field of view. With the 
advent of UKIDSS, UKIRT now operates in part 
as a survey telescope and in part as an open ac- 
cess telescope offering time through periodic peer- 
reviewed competition. At the time of writing, WF- 
CAM is scheduled for 60% of UKIRT time, 220 nights 
per year. After removal of engineering time, and time 
allocated to the University of Hawaii, and Japan, an 
average of 167 nights per year is left, 80% of which, 
134 nights per year, is devoted to UKIDSS, and the 
remaining time to other peer-reviewed programmes. 
(The latter are selected by the UK PATT system, 
open to world wide proposals). 

WFCAM is described in detail by Casali et al. 
(2007). Here we summarise some key characteris- 
tics. WFCAM has an unusual design, with an array 
of IR detectors inside a long tube mounted above 
Cassegrain focus. The forward-Cassegrain Schmidt- 
like camera design makes possible a very wide field 
of view (40 arcmin) on a telescope not originally 
designed for this purpose. The camera has four 
2048x2048 Rockwell Hawaii-II PACE arrays. The ar- 
rays have a projected pixel size of 0.4", which gives 
an instantaneous exposed field of view of 0.207 deg^ 
per exposure. The arrays are spaced by 0.94 detec- 
tor widths. The focal plane coverage is illustrated in 
Figure [5] 

The WFCAM filters, the optical performance, 
and the detector efficiency are presented in Casali et 
al. (2007). The photometric system is described in 
Hewett et al. (2006). 



3.2 Observing with WFCAM 

WFCAM makes short exposures on the sky, typically 
5-10 seconds. Coverage of the sky is then built up in 
several stages - by small jittering patterns around a 
fixed telescope pointing position; by macrostepping 
to make a filled in "tile"; and by accumulating sets 



* http://www.jach.hawaii.edu/UKIRT/ 
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Figure 5. The WFCAM focal plane. The spacing between 
detectors is 94% of the width of each detector, A sequence 
of four pointings therefore produces complete coverage for 
one "tile" , plus a small overlap region. The central dia- 
mond is the autoguider CCD. For details see Casali et al. 
(2007). 



of such tiles to gradually cover the sky, or revisiting 
tiles to build up depth. 

Integration overheads. Each WFCAM array has 
a separate SDSU controller, reading out 32 channels 
(8 segments in each of 4 quadrants). The total read- 
out time is 0.7 seconds, and the reset process is al- 
most instantaneous, but normal procedure involves 
two dummy readouts to stabilise the array, which take 
0.3 seconds. The overall exposure overhead is there- 
fore 1.0 seconds. Multiple exposures can be made at 
the same pointing position, co-adding into the same 
data file. The sum of these exposures is an integra- 
tion, which results in a set of four data files which in 
offline processing are later linked together as a multi- 
frame. For some users of WFCAM there may be sev- 
eral exposures per integration, but UKIDSS proce- 
dures always have only one exposure per integration. 
In addition each integration has a data acquisition 
overhead of 2-3 seconds. It is hoped that hardware 
and software improvements will improve this in due 
course, but for the moment there is then a total in- 
tegration overhead of 3-4 seconds. Exposures of 5-10 
seconds are therefore long enough to be reasonably 
efficient (~ 55-77%) and for sky background noise to 
be larger than the readout noise. The main exception 
is in the Y and Z bands, where exposures of 20 sec- 
onds are needed for the sky noise to exceed readout 
noise. Exposures longer than this are not normally 
used, as an increasing number of stars in the field are 
saturated. 

Jittering and micro-stepping. Small accurate 
telescope offset patterns relative to a fixed base po- 
sition are used to improve WFCAM data. The first 



method is to use a jitter sequence with offsets equal 
to whole numbers of pixels, resulting in frames which 
can be co-added. The aim of such a jitter pattern is 
to minimise the effects of bad pixels and other flat- 
fielding complications. A variety of jitter patterns can 
be used. The second kind of pattern is micro stepping, 
which uses offsets with non-integer numbers of pix- 
els. In 2x2 microstepping, offsets by N-l-1/2 pixels 
are used. The data are then interlaced (i.e. keeping 
the pixels independent) into a grid of pixel-spacing 
0.2", producing an image of size 4096x4096 pixels for 
each array. In 3x3 microstepping, offsets by N-l-1/3 
and N-f 2/3 pixels are used. The aim of such a mi- 
crostepping pattern is to improve image sampling - 
the WFCAM pixel size of 0.4" is adequate for mod- 
erate seeing (i.e. >0.8"), but undersamples the ex- 
pected seeing a significant fraction of the time. For 
sampling and/or cosmetic reasons, all UKIDSS sur- 
veys use at least 2 offset positions, and most use 4 
offset positions. The "standard shallow observation" 
has a total integration time of 40 seconds, with indi- 
vidual exposures chosen to add up to this total. The 
data from offset sequences are normally interleaved 
and/or stacked offline to make a single multi-frame 
data file, which is the basic unit of the archived data. 
Small offsets (less than 10") are intended to have min- 
imal overhead. They are applied by the tip-tilt sys- 
tem, with a smooth restoration to the optical axis 
during observing. Larger offsets are applied by the 
telescope mount, and observing commences after the 
system has stabilised. 

Tiling. The WFCAM arrays are spaced by 0.94 
detector widths. The sky could potentially be covered 
in a variety of mosaic patterns, but the typical proce- 
dure would be to expose in a pattern of four macro- 
steps to make a complete filled- in "tile" . Allowing for 
overlaps with adjacent tiles, the width of a single tile 
is then 3.88 detector widths, i.e. 0.883°, giving a solid 
angle of 0.78 sq.degs. The time between macrostep 
integrations (slew, stabilise, guide star lock) is ~ 15 
seconds. For shallow surveys, where each pointing has 
typically 4 offset positions each with a 10 second ex- 
posure and a 3.5 second overhead, this means that 
a tile with 4 pointings spends 160 seconds exposing 
out of an elapsed time of 276 seconds, making a to- 
tal observing efficiency of 58%. For deeper surveys, 
where many integrations are made between telescope 
slews, the macro-step overhead is negligible, and the 
efficiency tends towards ~75%. 

Schedule Blocks and Survey Definition. UKIRT 
operates an automated flexible queuing system. A 
precise sequence of exposures, offset patterns, and 
fllters at each of a list of pointing positions, which 
can be thought of as grouped into "tiles" as ap- 
propriate, is specifled in advance. These Obser- 
vations are grouped into "Minimum Schedulable 
Blocks (MSBs)", occupying roughly 20-60 minutes. 
The MSBs also contain constraints that determine 
whether they can be observed - required seeing, sky 
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brightness, etc. Calibration observations - twlight 
flats, standard stars, etc - are entered as indepen- 
dent MSBs. The MSBs are entered into a database 
which is queried during observing to generate a prior- 
ity ranked list of MSBs for which the current weather 
and observing conditions are suitable. The observer 
selects an MSB from this list (normally the highest 
priority MSB) and sends it onto an execution queue 
to be observed. If the MSB is sucessfuUy completed, 
it is marked as such in the database and will not be 
listed in future query results. (Occasionally MSBs are 
partially completed or fail the later Quality Control 
(QC) process, and so can be repeated). For UKIDSS, 
a Survey Definition Tool (SDT) is used to design the 
list of pointing positions and associated guide stars 
for each survey, which are then grouped into MSBs, 
and likewise into smaller "projects" which help plan- 
ning and monitoring of survey completion. 

Actual survey data collection rate in practice is 
considered in section [T] Implementation details, such 
as seeing and sky brightness constraints, are given 
in the data release papers - Dye et al. (2006) and 
Warren et al. (2007a,b). 

3.3 Survey calibration 

The UKIDSS data are calibrated to magnitudes in 
the Vega system. The WFCAM photometric system 
- filter response curves, and synthetic colours for a va- 
riety of objects - is described in Hewett et al. (2006). 
(Approximate passbands are given in section II. 2p . 
Calibration on the sky is achieved using observations 
of 2MASS stars within each field, which allows us 
to derive photometric calibration even during non- 
photometric conditions, including colour equations 
for transformation from the 2MASS system to the 
WFCAM system. For the shallow surveys, the Qual- 
ity Control (QC) process includes frames whose ze- 
ropoint is within 0.2 magnitudes of the modal value. 
(The vast majority of accepted frames are much bet- 
ter than this.) For the stacked surveys, where there is 
an additional self-calibration check, we accept frames 
within 0.3 magnitudes. It is possible that within such 
occasional frames with modest extinction there will 
be some noticeable spatial opacity gradient. This pos- 
sibility is not included in the current pipeline, but we 
intend to examine it in later re-processings. 

There are plenty of unsaturated 2MASS stars 
in every data frame - in the range 60-1000, de- 
pendendent on Galactic Latitude. Furthermore the 
2MASS global calibration is accurate to better than 
2% across the entire sky (Nikolaev et al. 2000). The 
procedure is to cross-match objects detected by the 
pipeline with 2MASS unsaturated sources that have 
(7jhk{'2MASS) ^ 0.1, and to transform the photom- 
etry of these stars into the WFCAM ZYJHK system 
using empirically derived colour terms. After correct- 
ing counts for the known radial variation in pixel 
scale, the average of these stars gives a global per- 



frame zeropoint. Tests against observations of UKIRT 
faint standards (Hawarden et al 2001) indicates that 
this procedure gives us a JHK photometric system ac- 
curate to 2%, which is the survey design requirement. 
(At the time of writing, the quality of the Z-, Y- and 
narrowband-filter calibration has not yet been quan- 
tified.) Calibration from 2MASS stars therefore seems 
justified. However, we have also made freqent obser- 
vations of UKIRT faint standards which provides a 
backup calibration, and an independent method of 
deriving colour equations. The calibration procedure, 
and the final colour equations between various sys- 
tems, will be presented in full in Hodgkin et al. (in 
preparation) . 



4 SURVEY DESIGN 

The design of the UKIDSS survey components was 
driven by a mixture of the legacy ambition, practi- 
cal limitations, and specific science goals. The total 
size of the project was chosen by a decision to con- 
tinue long enough to achieve a product of interna- 
tional significance and lasting value. The total time 
available was driven by UKIRT/WFCAM scheduhng 
constraints. We thus arrived at a seven year plan to- 
talling approximately a thousand nights. Seven years 
is however a long time to wait for science results; we 
therefore also designed an initial two year plan that 
would produce a self contained product and valuable 
science. 

Table [3] summarises the design parameters of 
each of the five UKIDSS survey components. Figure 
[6] shows the location of the survey fields on the sky. 
The logic behind this design, and some more detail 
about how the surveys are implemented, is described 
below for each survey component in turn. The im- 
plementation details were revised after the first ob- 
serving block (May- June 2005). We outline the cur- 
rent scheme, with the expectation that it is unlikely 
to change significantly. More complete details will be 
provided in each paper accompanying milestone data 
releases (Dye et al. (2006), for the EDR; Warren et 
al. (2007a,b) for DRl and DR2). Note that the sur- 
vey parameters shown here are the survey goals, as 
approved during the proposal stage, and adjusted to 
a minor degree during survey planning. Likewise the 
times quoted are the original proposal estimates. The 
actual time to reach these design parameters, or the 
likely final adjusted survey parameters, depend on 
both time allocation by the UKIRT Board and the 
real survey progress rate, which is discussed in sec- 
tion 



4.1 Microstepping strategy 

Microstepping improves the sampling, but has the 
disadvantage of making overheads worse. With 2x2 
microstepping, repeated after an offset, the exposure 
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0.77 
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Table 3. Summary of design goals for each survey. Note that the time estimates quoted are based on the efficiencies 
etc assumed at the time of proposal. The real time required to achieve these goals will be somewhat longer, and is being 
assessed at the end of the two year plan. Notes ; (i) Area is in square degrees, (ii) "Jx2" implies that two passes of 
the whole area are made in that filter. (Hi) "Limit" is the Vega magnitude of a point source predicted to be detected 
at 5(7 within a 2" aperture, (iv) tint *s the accumulated integration time at each sky position, (v) ttot is the number of 
hours required on-sky to complete the survey, allowing for the estimated exposure efficiency and mosaic efficiency. These 
efficiency factors are different for each survey, (vi) "Nights" is the estimated number of nights required, allowing for 
calibration and average fraction of useable UKIRT time (assumed 70% clear). 



time in the shallow surveys (mostly 40s total integra- 
tion) is 5s. Without microstepping longer exposure 
times are possible, 10s or 20s, which reduces the over- 
heads significantly. Experimentation with WFCAM 
data indicates that microstepping has little advan- 
tage for photometric accuracy, but improves the as- 
trometric accuracy. In each of the shallow surveys, 
repeat passes are therefore made in a particular 'as- 
trometric' band, for measuring proper motions; J for 
the LAS, and K for the GCS and GPS. As a compro- 
mise between increased overheads and better astrom- 
etry, in the current scheme for the high latitude shal- 
low surveys (LAS and GCS), the astrometric band 
is 2 X 2 microstepped, while the other bands are not 
microstepped. For the GPS, 2x2 microstepping is 
employed because of the importance of object sepa- 
ration in crowded fields. (Of course 3x3 would give 
even better PSF sampling, but this only makes a sig- 
nificant difference in the small percentage of cases 
where seeing is better than 0.3".) For the deep sur- 
veys, the total integration times in a field are much 
larger, allowing the use of longer exposures while mi- 
crostepping, so overheads are not an issue. In the DXS 
2x2 microstepping is used, and in the UDS 3x3 mi- 
crostepping is used, in all bands. 



4.2 Design of the Large Area Survey (LAS) 

The main science goals of the LAS require as large a 
volume as possible, with increased area being a more 
efficient use of time than increased depth. However 
the survey rate is constrained by our requirement 
for multiple pass bands, by the need for exposures 
long enough to avoid inefficient observing, and by the 
need for jittering in order to improve cosmetic quality 
and/or spatial sampling. 

The detailed implementation also depends on the 
weather constraints used. Because LAS uses a large 
fraction of the UKIDSS time, it would obviously not 
demand the best seeing. Colours are very important, 
and some objects are variable, which argues for doing 
all four bands when each sky position is visited; on the 
other hand, the Y and J observations require darker 
sky than the H and K observations. MSBs were then 
grouped so that H and K would be done together, and 
Y and J done separately. However, the queue is mon- 
itored and adjusted to try to make sure the YJ and 
HK observations are not too far apart. Our current 
plan is to prioritise uniformity and achieve the limits 
given in section (5] Thus in poorer seeing conditions 
the integration times are increased to compensate. 
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Figure 6. Location on the sky of the fields comprising the various survey components. Cross-hatch : Large Area Survey. 
Dark Grey : Galactic Plane Survey. Light Grey : Galactic Clusters Survey. Open rectangles : Deep Extragalactic Survey. 
Note that the Ultra Deep Survey lies just to the west of the DXS field at 02hl 8m - 05° 1 ' The dashed line marks the 
Galactic plane. Note that UKIRT lies at latitude +20° . 



Field selection for the LAS was designed to have 
a good spread in RA, to have a reasonable amount 
of sky coverage at lower declinations, for follow-up 
on ESO telescopes, to keep below the UKIRT Dec- 
lination limit (-f60°), all while lying with the SDSS 
footprint. There are three sub-areas, shown in Fig 
[B] The detailed field co-ordinates are defined on the 
UKIDSS web pages. 

{i)The LAS equatorial block : 1908 sq.deg. This 
includes most of SDSS stripes 9 to 16. 

(u) The LAS northern block : 1908 sq.deg This 
includes most of SDSS stripes 26 to 33. 

(iii) The LAS southern stripe : 212 sq.deg. This 
is a section of SDSS stripe 82, extending over - 
25° <RA<-h60°, -1.25° <Dec<-f 1.25°. Stripe 82 has 
been repeatedly scanned by SDSS, and this is the re- 
gion of highest quality. 



4.3 Design of the Galactic Plane Survey 

(GPS) 

The GPS aims to map as much of the Galactic Plane 
as possible to a latitude of ±5°. The Galactic Lati- 
tude limit is chosen to match other surveys, for ex- 
ample the MSX survey (Egan and Price 1996). The 
survey area is then largely dictated by the UKIRT 
Declination limit to the North, and by the latitude 
of UKIRT to the South. Within a reasonable length 
of time, we can then afford to go roughly a factor of 
two deeper than the "standard shallow observation" 
defined in section 13.21 This depth is good enough to 
see all of the IMF to the H-burning limit in quite dis- 



tant clusters, to see AGB stars all the way through 
the Galaxy, and to see ordinary G-M stars to several 
kpc. The tradeoff to consider is then between depth, 
colours, and repeat coverage. At least three bands, 
and preferably, four are needed, in order to estimate 
both spectral type and extinction. However, extinc- 
tion is large enough in much of the Plane that Y band 
observations are impractical. For regions of low ex- 
tinction the optical IPHAS survey at r',i' and Halpha 



(see http://www.iphas.org I will provide sufficient ad 



ditional colours to determine the average extinction 
as a function of distance in each field using reddening 
independent colour indices. Other statistical methods 
to measure extinction using the JHK colours alone 
can also be employed - see Lopez-Corredoira et al. 
(2002). Measurement of both variability and proper 
motions is a goal of the GPS. As a bare minimum 
to achieve this, we plan three epochs in one band 
spread across ~ 5 years. We choose the K-band to 
make these repeats, because, given extinction, this is 
the sensible band in which to build up depth - it is 
the K-band that allows us to see clean through the 
Galaxy in most directions. In summary then, we plan 
an initial pass at JHK, with integrations of 80s, 80s 
and 40s in the 3 bands respectively, followed by two 
further passes at K with 40s integrations at intervals 
of at least 2 years for any survey tile. 

One of the goals of GPS is the discovery and 
study of Young Stellar Objects, and in particular 
molecular outflows. We therefore plan in addition to 
the above a survey of a single large star forming re- 
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gion in a narrow band H2 filter. (For details of tliis 
filter, see Casali et al. 2007). 

The area to be mapped is shown in Fig. [H] The 
main area is defined by the Galactic latitude range 
b = ±5°, Dec<60°, and Dec >-15°. These constraints 
define two sections of Galactic longitude, which are 
15° < i < 107°, and 142° < I < 230°. In addi- 
tion we will map a narrow extension through the 
Galactic centre, within b = ±2°, covering Galac- 
tic latitudes -2° < / < 15°. The Galactic bulge 
will also be explored by surveying a thin stripe ex- 
tending upwards in latitude from the Galactic cen- 
tre. Finally, the molecular hydrogen survey maps the 
Taurus- Auriga-Perseus complex. 

Normal procedure is to do all bands in one visit, 
as colours are important and many objects are vari- 
able. 



4.4 Design of the Galactic Clusters Survey 
(GCS) 

The design of the GCS is relatively simple. It needs to 
target several separate clusters, in order to examine 
the substellar IMF over a range of ages and metal- 
licities. The depth requirement is set by the need to 
detect objects in the 30-100 Mj range. Assuming the 
standard "shallow survey" depth (see section I3.2|l . 
this means that clusters have to be fairly close and/or 
young, i.e. within a few hundred parsecs and/or less 
than a few hundred million years old. There are rela- 
tively few such objects, and they are several degrees 
across. A natural strategy therefore emerges using 
the standard shallow depth and surveying ten nearby 
clusters, covering 1067 sq. deg. in total. To distin- 
guish cluster members, all five passbands are needed 
(ZYJHK), plus a measurement of a proper motion 
using a second pass in the K band. Full colour in- 
formation provides cluster sequence discrimination in 
multi-colour space, reddening estimates using shorter 
versus longer wavelength colour indices, breaking the 
degeneracy between reddening due to instellar extinc- 
tion and that due to the presence of circumstellar 
disks. 

The strategy is to cover the majority 1067 sq. 
degs in a single pass in K, to provide the proper mo- 
tion baseline, within the initial two year plan. Table 
|4] lists the parameters of the chosen clusters. 

4.5 Design of the Deep Extragalactic 
Survey (DXS) 

The DXS aims to detect galaxies at redshifts of 1-1.5. 
To avoid selecting only the brightest and hence most 
massive galaxies, this requires the detection of galax- 
ies close to the break in the galaxy luminosity func- 
tion, Mk which is -22.6 locally (Bell et al. 2003). At 
z=l this corresponds to an evolution corrected, total 
K magnitude of 20.7, and 21.8 at z=1.5. Therefore, 
taking into account aperture effects, our target depth 



of K=21 will reach to within 0.5-0.7 magnitudes of 
Mk and hence sample a representative galaxy popu- 
lation at z=l. The NIR galaxy colours at this redshift 
lie in the range J-K=1.5-1.8 so to provide a photo- 
metric constraint on the galaxy redshift we also re- 
quire observations to J=22.3 to ensure matched J and 
K detections for the target galaxies. 

The survey area was driven by the aim to sam- 
ple large scale structure at z=l on scales and volume 
comparable to that measured locally (~100 Mpc and 
0.2 Gpc'' respectively). At z=l, our assumed cosmol- 
ogy implies that 100 Mpc corresponds to 3.5° and 
a 0.2 Gpc^ volume in the range z=l-1.5 requires 
40 sq.deg. Therefore a minimum combination of 3x3 
WFCAM tiles will span these scales and a total of 54 
WFCAM tiles would be required to cover that area 
(0.75 sq.deg. per field). 

The number and position of the DXS survey 
fields were chosen to provide the best combination of 
quality and coverage of supporting, multi- wavelength 
data, to maximise the spatial scale sampled by each 
individual field (~ 100 Mpc) for clustering studies 
and allow a uniform coverage in right ascension. Bal- 
ancing these factors resulted in the selection of four 
survey fields: 1) XMM-LSS (centre: 2h25m -04d30m) 
- a SWIRE, CFHTLS, VVDS, GALEX and XMM 
survey field adjacent to the UKIDSS UDS area; 2) the 
Lockman Hole (centre: 10h57m -|-57d40m) - centred 
on the SHADES survey area but within the SWIRE 
and GALEX survey areas with extensive radio cov- 
erage; 3) ELIAS-Nl (centre: 16hl0m +54d00m) - 
a SWIRE and GALEX field with additional radio, 
optical and X-ray data; 4) SA22 (centre: 22hl7m 
+00d20m) - centred on VVDS-4 but the least well 
surveyed area, included to ensure a uniform demand 
with RA. Other survey fields were considered (COS- 
MOS, NOAO-DWFS, Groth Strip, Spitzer-FLS) but 
most were either not sufficiently large or comprehen- 
sive to justify inclusion. The chosen fields are listed 
in Table [5] 

The total area covered by the DXS within the 
full 7 year span of UKIDSS will depend on weather 
and competition from other surveys (most notably 
the UDS) but our goal is 35 sq.deg. or 12 WFCAM 
fields in each DXS survey area in J and K. We also 
intend towards the end of the survey to include an 
additional 1-2 WFCAM fields in the centre of each 
DXS survey area in H to broaden the photometric 
coverage. 

For DXS, star-galaxy separation at faint magni- 
tudes will be very important, so 2x2 micro-stepping 
is employed to give good sampling. Achieving the 
required depth will require reliable stacking, and so 
minimising any systematic effects in detector struc- 
ture that do not fiat-field out. The DXS strategy 
therefore employs substantial jittering. Each visit to 
a given tile position uses 10 second exposures, a six- 
teen point jitter, and 2x2 microstepping at each of 
these jitter positions. Each such visit therefore has 
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Priority/ 


Type 


RA 


Dec 


Area 


Age 


Minimum 


Name 




(2000) 


(2000) 


(sq.deg.) 


(Myr) 


mass (Mq) 


(1) IC 4665 


open cl. 


17 46 


+05 43 


3.1 


40 


0.020 


(2) Pleiades 


open cl. 


03 47 


+24 07 


79 


100 


0.024 


(3) Alpha Per 


open cl. 


03 22 


+48 37 


50 


90 


0.025 


(4) Praesepe 


open cl. 


08 40 


+19 40 


28 


400 


0.046 


(5) Taurus-Auriga 


SF assoc. 


04 30 


+25 00 


218 


1 


0.010 


(6) Orion 


SF assoc. 


05 29 


-02 36 


154 


1 


0.014 


(7) See 


SF assoc. 


16 10 


-23 00 


154 


5 


0.010 


(8) Per-OB2 


SF assoc. 


03 45 


+32 17 


12.6 


1 


0.011 


(9) Hyades 


open cl. 


04 27 


+15 52 


291 


600 


0.041 


(10) Coma-Ber 


open cl. 


12 25 


+26 06 


79 


500 


0.043 



Table 4. Areas targeted for the Galactic Clusters Survey, listed in priority order. The primary data source is Allen (1973) 
with updated Hipparcos distances from Robichon et al (1999) and de Zeeuw et al (1999); assuming a depth J=19.5 (Table 
[M the mass limit was calculated using the DUSTY models of Chabrier et al (2000). 



an integration of 640 second at each sky point. To 
reach the intended depth requires a total exposure of 
2.1 and 1.5 hours in J and K respectively or 12 and 
8 visits each. H observations are quite slow, and are 
planned as a lower priority in a subset of the chosen 
fields. Given that the DXS observations do not re- 
quire photometric conditions or the very best seeing, 
then the final number of visits for each field may be 
higher to compensate for these poorer conditions. 



4.6 Design of the Ultra Deep Survey (UDS) 

The UDS aims to go as deep as possible in a sin- 
gle contiguous WFCAM tile. The depth is set by the 
aim of detecting giant ellipticals at z—3 if they exist. 
The total magnitude of such objects is expected to 
be K~ 21 but they will be significantly extended, so 
that we need to reach a point source depth of K=23. 
Three bands are needed, to get photometric redshifts 
and discriminate between objects. To effectively sep- 
arate ellipticals and starbursts, we need to be able to 
detect colours J-K~ 2 and H-K~l, otherwise most of 
our detections may be K-band only. This sets limits 
of J~25 and H~24, which are in fact more demand- 
ing in time than the K-band observations. The final 
expected depths are J=24.8, H=23.8, K=22.8. 

As with the DXS, we need good sampling to 
enable star-galaxy separation at faint magnitudes, 
and multiple jitters to overcome detector systematics 
when stacking. These issues are even more demand- 
ing however for the UDS; at each visit we use 3x3 
microstepping and a 9 point jitter, and repeated vis- 
its are not at precisely the same position, but in a 
carefully arranged pattern - a kind of super-jitter. 

Although the headline science goals of the UDS 
are to find high redshift galaxies which will be ex- 
tremely faint at optical wavelengths, the legacy value 
of the survey will be increased by the presence of com- 
plementary ultra-deep optical data. The field chosen 
was therefore the Subaru/XMM-Newton Deep Field 
(Sekiguchi et al. , in preparation) which possessed the 



deepest optical data of any square-degree field at the 
time this decision was made. Considerable data at 
other wavelengths from radio to X-ray also exist in 
this field. This field is at the western edge of one of 
the DXS fields. 



4.7 Two Year Goals 

We aim to complete self contained and scientifically 
valuable datasets on a two year timescale. The de- 
tailed plan is set out in Dye et al. (2006), but briefly is 
as follows. The shallow surveys (LAS, GCS, and GPS) 
are accelerated compared to the deep stacked surveys 
(DXS and UDS). In addition, the LAS concentrates 
on southern latitudes in the first two years, in order to 
maximise VLT follow-up. For LAS, roughly half the 
area - the equatorial block, and the southern stripe 
- will get complete YJHK coverage. (Second epoch 
J for the same areas will come later). Additional J 
only coverage in the Northern block will be achieved 
as time permits. For GPS, the prime aim is to obtain 
the first of three K epochs over the whole survey area, 
with J, H, and H2 coverage over a sub-area. For GCS 
the aim is to get complete filter coverage for five of the 
ten target clusters, and the central regions of three, 
plus K only coverage of the remaining two. For DXS, 
the two year aim is to reach the full depth in J and 
K for a subset of the area - four tiles (3.1 square de- 
grees) in each of the four fields. The UDS is of course 
a single tile. The two year goal is to achieve K=22.8 
(full depth) and J=23.8 (one magnitude short), with 
no H coverage. 



5 DATA PROCESSING AND DATA 
PRODUCTS 

The commitment to making a public survey requires 
the construction of complete, reliable, tested, and 
documented products from the raw data. The very 
large volume of UKIDSS data (200 GBytes/night) 
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Survey Area RA(2000) Dec (2000) 



XMM-Subaru 


UDS 


0.77 


02 18 00 


-05 10 00 


XMM-LSS 


DXS 


8.75 


02 25 00 


-04 30 00 


Lockman Hole 


DXS 


8.75 


10 57 00 


-1-57 40 00 


ELAIS Nl 


DXS 


8.75 


16 10 00 


-1-54 00 00 


SA22 


DXS 


8.75 


22 17 00 


-f 00 20 00 



Table 5. Fields targeted for the Deep Extragalactic Survey and the Ultra Deep Survey. Note that the UDS field is at the 
western edge of the DXS XMM-LSS field 



means that to achieve these goals requires a uni- 
form and automated approach to data processing. 
Likewise the large accumulated volume of prod- 
ucts - expected to be several tens of Terabytes of 
image data and several billion source detections - 
means that as well as providing public data ac- 
cess, we need to provide online querying and anal- 
ysis facilities as a service. These ambitions are met 
for all WFCAM data (both UKIDSS survey data 
and PATT PI data) by the VISTA Data Flow Sys- 
tem (VDFS). VDFS is a PPARC- funded project in- 
volving QMUL, Cambridge and Edinburgh, aimed 
at handling the data from first WFCAM and then 
the VISTA telescope. (The Science Archive was also 
prototyped on the SuperCosmos Science Archive : 
see http://surveys.roe.ac.uk/ssa/ and Hambly et al. 
2004). The system aims at (i) removing instrumen- 
tal signature; (ii) extracting source catalogues on a 
frame by frame basis; (iii) constructing survey level 
products - stacked pixel mosaics and merged cata- 
logues; (iv) providing users with both data access and 
methods for querying and analysing WFCAM data. 
Overall data flow is as follows. Raw data are 
shipped by tape on a weekly basis from Hawaii to 
Cambridge, where they are available within a month 
of the observations being taken. Raw data are then 
transferred via the internet for ingest into the ESO 
archive system. Pipeline processed single frame data 
are transferred to Edinburgh over the internet on a 
daily basis, where they are ingested into the science 
archive, and further processing (stacking, merging, 
and quality control) takes place. The processed data 
are then released to the public at periodic intervals. 



5.1 The WFCAM Pipeline 

The general philosophy behind the pipeline process- 
ing is that all fundamental data products are FITS 
multiextension files with headers describing the data 
taking protocols in sufficient detail to trigger the ap- 
propriate pipeline processing components, and that 
all derived information, quality control measures, 
photometric and astrometric calibration and process- 
ing details, are also incorporated within the FITS 
headers. Generated object catalogues are stored as 
multiextension FITS binary tables. These FITS files 
thereby provide the basis for ingest into databases 
both for archiving and for real time monitoring of 
survey progress and hence survey planning. 



After conversion at the summit from Starlink 
NDF to FITS files, to reduce the data storage, I/O 
overheads and transport requirements, we make use 
of lossless Rice tile compression (eg. Sabbey et al. 
1998). For this type of data (32 bit integer) the Rice 
compression algorithm typically gives an overall fac- 
tor of 3-4 reduction in file size. Data are shipped 
roughly weekly from JAC using LTO tapes, one per 
detector channel, and combined to create the raw 
archived multiextension FITS files on ingest in Cam- 
bridge. 

The data processing strategy attempts to min- 
imise the use of on-sky science data to form "calibra- 
tion" images for removing the instrumental signature. 
By doing this we also minimise the creation of data- 
related artefacts introduced in the image processing 
phase. To achieve this we have made extensive use of 
twilight fiats, rather than dark-sky flats (which po- 
tentially can be corrupted by thermal glow, fringing, 
large objects and so on) and by attempting to decou- 
ple, insofar as is possible, sky estimation/correction 
from the science images. 

Each night of data is pipeline processed indepen- 
dently using the master calibration twilight flats (up- 
dated at least monthly) and a series of nightly gen- 
erated dark frames covering the range of exposure 
times and readout modes used during that night. A 
running sky "average", i.e. calculated from a number 
of recent frames, in each passband is used for sky arte- 
fact correction. After removing the basic instrumental 
signature the pipeline then uses the header control 
keywords to produce interleaved and/or combined 
(stacked) image frames for further analysis. This in- 
cludes generation of detected object catalogues, and 
astrometric and photometric calibration based on 
2MASS. A more detailed description of the WFCAM 
processing is given in Irwin et al. (in preparation). 



5.2 The WFCAM Science Archive 

Data processing delivers standard nightly pipeline 
processed images and associated single passband cat- 
alogues, complete with astrometric and flrst-pass 
photometric calibrations and all associated 'meta' 
(descriptive) data in flat FITS files. These data are 
ingested into the archive on a more or less daily basis. 
To produce UKIDSS survey products however three 
more processes are needed - image stacking, source 
merging, and Quality Control (QC) filtering. 
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Source merging is computationally very inten- 
sive, UKIDSS catalogues are extremely large, and the 
VDFS system is designed to scale to the even larger 
expected VISTA products. For standard UKIDSS 
products therefore a relatively simple spatial pair- 
ing algorithm is used. A pairing radius of roughly 
1" (varying between sub-surveys) is used, several 
times larger than the astrometric accuracy (0.1"), 
in order to allow merging of moving objects between 
epochs. Note that the WSA itself foUows SDSS prac- 
tice in storing a neighbour table for every source, so 
that much more flexible source matching algorithms 
can be applied later, which may be crucial when fol- 
lowing up objects with rare colours and so on. More 
details of stacking, mosaicing, and source merging are 
given in Irwin et al. (in preparation) and Hambly et 
al. (in preparation). 

The QC process is a joint responsibility of the 
UKIDSS consortium and the VDFS project and is a 
semi-automated process. The first aim is to remove 
data that is corrupted or unusable (for example due 
to bright moon ghosts) . This is achieved partly by au- 
tomated inspection of headers and data, and partly 
by visual inspection. The second aim is to deprecate 
a subset of data frames not achieving a standardised 
'survey quality' by applying a series of cuts based on 
QC parameters derived from the images - for exam- 
ple, seeing, average stellar ellipticity, sky brightness 
and so on. These cuts are evolving with experience at 
each data release. Histograms of QC parameters are 
derived, and manual decisions taken on the cut lev- 
els, which are then applied automatically. Each sub- 
survey (LAS, UDS etc) has separate QC parameters. 
For example, in the EDR, the LAS seeing limit was 
1.2", whereas for the UDS it was 0.9"and 0.8"at J and 
K respectively. The process is described more fully in 
the "Early Data Release (EDR)" paper of Dye et al. 
(2006), and revised parameters for DRl are provided 
in Warren et al (2007a). 

Image data volume is typically ~200 GBytes per 
night, with catalogue and descriptive data being typ- 
ically ~ 10% of that figure. Hence, over the course 
of several years of observations it is anticipated that 
10s/ 100s of Tbytes of catalogue/image data will be 
produced by survey operations with WFCAM. In or- 
der to enable science exploitation of these datasets, 
the concept of a 'science archive' has been developed 
as the final stage in the systems-engineered data flow 
system from instrument to end-user (Hambly et al. 
2004b). 

Finally, we note that the WSA has been devel- 
oped throughout with Virtual Observatory compat- 
ibility in mind. The WSA exports VOTABLE, ex- 
poses its query interface as an ADQL web service, 
is published in VO registries, and provides access to 
AstroGrid VO tools such TopCat and MySpace. 
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Figure 7. The Y-J,J-H two colour diagram for a sin- 
gle tile observed in the LAS SV programme. Black dots 
show the data for stellar sources detected in the WFCAM 
data. Also shown are the synthetic colours of QSOs, L/T 
dwarfs, and model Y dwarfs. The solid line shows the po- 
sitions of M dwarfs. The observation was targeted at a 
known T2 dwarf SDSS J125454-012247 which is recov- 
ered with Y-J=1.10 and J-H=0.54. 



The WFCAM Science ArchiveQ (WSA) is much 
more than a simple repository of the basic data prod- 
ucts described previously. A commercial relational 
database management system (RDBMS) deployed on 
a mid-range, scalable hardware platform is used as 
the online storage into which all catalogue and meta 
data are ingested. This RDBMS acts as the backing 
store for a set of curation applications that produce 
enhanced database driven data products (both image 
products, e.g. broad-band/narrow-band difference 
images; and catalogue products, e.g. merged multi- 
colour, multi-epoch source lists). Moreover, the same 
relational data model is exposed to the users through 
a set of web-interface applications that provide ex- 
tremely flexible user access to the enhanced database 
driven data products via a Structured Query Lan- 
guage interface. The primary purpose of the WSA 
is to provide user access to UKIDSS datasets, and a 
full description, along with typical usage examples, 
is given in Hambly et al. (in preparation). Step-by- 
step examples of WSA usage are also included in the 
UKIDSS EDR paper. Dye et al. (2006). 



^ http://surveys.roe.ac.uk/wsa 
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Figure 8. Illustration of the KX method from some 
10 sq degs of science verification observations in high 
Galactic latitude fields. All detected objects in the range 
14-5<K<16.5 are plotted, totaling 21000 sources. Stars 
make up the long, thin cloud, and galaxies form the 
shorter cloud to the right. The solid squares are model 
quasar colours 0<z<8, Az=0.1, from Hewett et al. 
(2006), with z=0 marked by the open square. Similarly 
triangles mark the model colours of an unevolving ellip- 
tical galaxy 0<z<3, Az=0.1. The large filled circles are 
the 11 SDSS spectroscopically confirmed quasars in the 
observed fields, brighter than the fainter limit of K=17. 
The diagonal dashed line represents a possible KX selec- 
tion criterion. Candidate quasars (of which there are 167) 
are compact sources to the right of the line, and are in- 
dicated by the open circles. Reddening vectors at different 
redshifts run approximately parallel to this line (Warren 
et al. 2000). 




Figure 9. Central 25% of a GPS K-band tile pointed at 
M17, showing the richness of data in the GPS. The full 
tile contains 740.305 sources. 



6 SCIENCE VERIFICATION 
PROGRAMME 

Because the consortium has no proprietary rights, 
and no formal role in deriving science results from the 
data (as opposed to facilitating their exploitation) , in 
this paper we have avoided showing science results 
from the released data. Examples of early papers us- 
ing UKIDSS data are given in section 17.31 However, 
the potential of UKIDSS can be well illustrated using 
the Science Verification data. 

Following the technical commissioning of WF- 
CAM, and before the commencement of formal survey 
operations in May 2005, the UKIDSS consortium un- 
dertook a modest set of test observations, as a "Sci- 
ence Verification (SV)" programme. These observa- 
tions were aimed primarily at further technical com- 
missioning, testing and tuning the implementation 
strategy, and exercising the data fiow system. How- 
ever the data collected have clearly demonstrated the 
scientific power of UKIDSS, and confirm the efficacy 
of the survey design. In this section we show some 
examples of science results from these SV data. 



6.1 Science Verification results for the LAS 

The LAS SV programme covered some 20 square de- 
grees, achieving close (0.2mag) to the standard shal- 
low depth in filters Y,J,H,K. One scientific aim was 
to test the likely recovery of cool brown dwarfs. The 
success of this is illustrated in Fig. [T] which shows 
data from a tile aimed at a known T2 dwarf, SDSS 
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J125454-012247 (Knapp et al. 2004). This object was 
indeed detected, and the colours found are consis- 
tent with those previously published. Other L and T 
dwarfs were also targeted and successfully detected. 
In this very limited area no new objects of significant 
interest have been found but that is in line with ex- 
pectations. The LAS SV data were also used to verify 
the photometry and astrometry of the UKIDSS sur- 
vey. Details of these tests are given in the EDR data 
release paper, Dye et. al. (2006). 

A second aim was to test the location of quasars 
by combining UKIDSS and SDSS colours. Figure [5] 
shows that the "K excess" method works extremely 
well. The stellar and galaxy sequences are cleanly 
distinguished. Point-like objects to the right of the 
dashed line are good quasar candidates. Known SDSS 
quasars in these fields are in the upper part of this 
region, but the UKIDSS SV data shows many more 
candidate quasars with similar colours. Addition- 
ally there are several candidates with much redder 
colours than the known quasars. These are candi- 
date reddened quasars, and spectroscopy is required 
to investigate their nature. Several have colours simi- 
lar to galaxies, but the overall colour spread of the 
candidates is much broader than for the galaxies. 
The quasars in these fields, if confirmed, will be at 
relatively modest redshift. Note that quasars move 
rapidly redder in g-i beyond z~3.8. The very high 
redshift quasars that we hope to find will be much 
sparser on the sky, and in gJK will be hard to distin- 
guish from cool brown dwarfs. Here, as explained in 
section [231 the Y-J colour will be crucial. Analysis of 
relevant data is still in progress and will be reported 
in a later paper. 

6.2 Science Verification results for the GPS 

The science verification data for the Galactic Plane 
Survey included a 0.75 deg^ tile centred on M17 (a 
high mass star formation region). The Ml 7 region 
has been well studied (eg. Jiang et al. 2002) and it 
provided a good test of the photometric reliability of 
the data in a nebulous region. 

The central 25% of the M17 tile is shown in Fig- 
ure |9] This illustrates the sensitivity of WFCAM to 
the structure and stellar population of distant star 
formation regions. By visual comparison of the WSA 
catalogue and the reduced image, it is clear that very 
few spurious sources are found in the archive but 
within the brightest nebulosity the archival source 
lists become seriously incomplete, with almost all 
point sources being undetected. This problem occurs 
generally throughout the GPS in very bright nebulae, 
since the UKIDSS pipeline photometry is designed to 
detect resolved galaxies as well as point sources, and 
hence it attempts to interpret bright nebulosity as 
a single extended source. Such very bright nebulae 
represent a tiny fraction of the total survey area. In 
general, more complete luminosity functions may be 



derived in nebulous regions by performing indepen- 
dent photometry on the reduced data, using the zero 
points provided in the image headers. 

At the present time only aperture photometry is 
available in the WFCAM Science Archive. One might 
expect the photometric precision to suffer consider- 
ably from the effects of source confusion. However 
the algorithm simultaneously fits the coordinates and 
fluxes of adjacent stars in small groups and divides up 
the fluxes between any overlapping apertures in a sen- 
sible manner (following principles described by Irwin 
1985). Hence the resulting photometric quality is very 
good even in very crowded flelds with source densities 
in excess of 10® per square degree, and is similar to 
the quality of results from proflle fitting photometry, 
which is a far more time consuming process. (Detailed 
comparison of various photometry techniques will be 
given in Irwin et al. (in preparation)). 

The quality of the photometry is illustrated in 
the two colour-colour diagrams in Figure [TOl each of 
which shows those sources within a single detector 
array (13 arcminute field) that are listed as having 
well measured photometry (errors < 0.05 mag on each 
axis) in the WSA. The curves shown are for unred- 
dened luminosity class V main sequence stars (lower 
curve) and luminosity class HI giants (upper curve), 
using synthetic photometry from Hewett et al. (2006). 
The left hand panel shows a field at (/, fe)=171.2, 4.8, 
where there is very little extinction. Nearly all of the 
stars lie near the curves, and both the G-type to 
early M type main sequence and the giant branch 
are well represented. There are relatively few sources 
with colours consistent with those of late M dwarfs 
near the end of the class V sequence. Hence bona 
fide late M and L dwarfs should be detectable by fol- 
low up observations with a reasonable success rate, 
provided they have precise photometry. The task will 
become easier when proper motion information be- 
comes available with the 2nd and 3rd epoch data. 

The right hand panel of [TU] illustrates a crowded 
field with high extinction at (Z,6) = 15.0, 0.0. Despite 
significant source confusion, there is a clear sepa- 
ration between the reddened giant branch and the 
reddened main sequence. The majority of the stars 
with low reddening are K-type to early M type dwarf 
stars. Giant stars are generally more distant and suf- 
fer higher reddening, so the giant branch becomes 
clearly defined only at {J - H) > 1.2, {H - K) > 0.6. 
Some curvature of the reddening sequences is appar- 
ent, which we attribute to the change in the effec- 
tive wavelengths of the JHK bandpasses at high red- 
dening. The number of photometric outliers can be 
greatly reduced by using additional quality informa- 
tion in the archive, eg. using source ellipticity to re- 
move binaries which are unresolved in one or two of 
the passbands. This will be explored fully in a later 
paper. 

Since the main sequence curve is almost parallel 
to the reddening vector it appears that additional col- 



© 2007 RAS, MNRAS 379, 1599-1617 



The UKIRT Infrared Deep Sky Survey (UKIDSS) 19 



1.5 



0.5 




- 


1 ' ' 








- 


Heddfciiint. 


Line 




- 


- 


/ 




•f 


- 


- 




< 




- 


r 




r 




~ 


_ 




f 




_ 








, 


1 1 1 



H-K 



H-K 



Figure 10. Two colour diagrams for two individual array detectors (12.8 arcmin fields), showing sources listed as having 
errors <0.05 mag on each axis in the WSA. Left hand panel ; an uncrowded field with low extinction at (l,b)=171.2,4-8, 
showing 1236 sources. Right Hand panel ; a crowded field with high extinction at (l,b)=15.0,0.0, showing 12954 sources. 
In each panel the lower curves shows the locus of the main sequence, and the upper curve shows the locus of luminosity 
class III giants. Even in the crowded field the photometric precision is good enough to show the 0.25 mag vertical separation 
between the reddened giant sequence and the reddened main sequence stars. An approximate reddening vector is plotted 
in the right hand panel. Some curvature of the reddening sequences is seen in the data, which is attributed to the change 
in the effective bandpasses at high extinction. 



ors will be required to permit precise photometric de- 
terminations of source extinction and hence spectral 
type and luminosity class. For blue sources this will 
be done with the aid of optical data from the IPHAS 
survey (www.iphas.org) while for very red sources 
mid infrared data from the SPITZER-GLIMPSE sur- 
vey (www.astro.wisc.edu/sirtf/) of part of the galac- 
tic plane and later the NASA WISE survey of the 
whole sky (wise.ssl.berkeley.edu/news.html.) may be 
useful. 

Profile fitting photometry will be included in a 
future release of UKIDSS data if it is found to im- 
prove the photometric precision. Experiments with 
DAOPHOT in IRAF and a preliminary profile fitting 
scheme for the microstepped data (D. Wyn Evans, 
private comm.) indicate that it is unlikely to produce 
a major improvement over aperture photometry but 
it may be possible to reduce the number of photomet- 
ric outliers and thereby aid the detection of sources 
with unusual colours. 

In relatively uncrowded regions the aperture pho- 
tometry reaches close to the target depths at J, H 
and K (see results in Warren et al. 2007). Source 
confusion does not significantly affect the measured 
depth (K=18.1) in the outer galaxy section of the 
GPS at ^=142-230. Even in the inner galaxy section 
(—2 < I < 107) the reduction in sensitivity is <0.5 
mag in most fields. However the sensitivity suffers 
more in the most crowded parts of the mid-plane as 
the Galactic Centre is approached, eg. it is reduced 
by 1.5 mag at 1=15, b=0. 



A full analysis of the data quality in the GPS 
and some results from the SV data will be presented 
in a future paper (Lucas et al. in preparation). 



6.3 Science Verification results for the GCS 

Science verification observations for the GCS yielded 
8 tiles in each of the three targets observable at that 
time: IC 4665, Upper Scorpius and Coma Berenices. 
In the case of IC 4665, the SV observations complete 
the required survey for that cluster. In Figure \TT\ we 
show a Z versus Z-J colour-magnitude diagram for 
Upper Scorpius. The observations cover 6 square de- 
grees and approximately 100,000 point sources are de- 
tected. Nearly all of these are of course background 
stars, with perhaps some foreground stars. In Fig llll 
we show only those sources with J>10.5 and Z>11.5. 
The main sequence and giant branch are clearly seen. 
In addition one can notice a clean sequence to the 
right of the diagram running from (Z-J, Z) = (l. 0,12.0) 
to (Z-J,Z) = (2.5,19.0), which must be the cluster se- 
quence. Candidate cluster members, plotted as large 
points in Fig. 1111 have been isolated first by selection 
in the Z-J vs Z diagram, and then in Y-J vs J-K to 
eliminate field dwarfs. Objects brighter than J=15.8 
are also detected in 2MASS, giving a proper motion 
estimate which allows a very reliable final determi- 
nation of cluster membership. For fainter objects we 
have only the photometric information, but estimate 
that field contamination is less than 10%. (See Lodieu 
et al. (2007) for details). 
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Figure 11. Z vs Z-J diagram for six square degrees in the 
Upper Sco field, showing only those sources with J>10.5 
and Z>11.5. Candidate cluster members have been iso- 
lated by colour cuts in multicolour space (see text). One 
in ten of the field sources are plotted, and all the clus- 
ter member candidates. Model isochrones are described in 
the text. The cluster sequence can be clearly seen following 
the models well. Objects on or to the right of the model 
isochrones, and so likely to be cluster members, are plotted 
as large points. The expected positions of objects of var- 
ious masses are indicated. Sub-stellar objects are clearly 
detected. 



Overplotted are 5 Myr theoretical isochrones 
shifted to a distance of 145pc, appropriate for the 
estimated age and distance of Upper Sco : BCAH98 
or NextGen models (solid line; Baraffe et al. 1998), 
DUSTY or BCAHOO (dashed line; Chabrier et al. 
2000), and COND03 (dotted line; Baraffe et al. 2003). 
These isochrones were specifically computed for the 
WFCAM fihers (Isabelle Baraffe and France AUard; 
personal communication). On these isochrones we 
also indicate object masses. We can clearly see there- 
fore that we are indeed locating brown dwarfs within 
the cluster and will be able to derive a substellar mass 
function all the way down to 10 Jupiter masses. 



6.4 Science Verification results for the DXS 
and UDS 

The DXS and UDS undertook a joint science veri- 
fication programme to establish the performance of 
WFCAM for deeper exposures. Fig 1121 shows the K- 
band number counts from an accumulated exposure 
of 1.5 hours (6 hours of data) on a single tile in the 
ELAIS-Nl field. Simulations show that the 50% com- 
pleteness limit in this field is K=20.3, significantly 
shallower than the expected 5a limit in this time. 



This is partly but not wholly because of poorer than 
average seeing in the SV observations. The reason for 
this discrepancy is discussed in section [7] Note that 
this figure includes both stars and galaxies. For com- 
parison, we show number counts from several other 
surveys. This figure illustrates the power of UKIDSS, 
as we can determine number counts accurately from 
a single tile. When the survey is complete, we will 
therefore be relatively immune to cosmic variance. 

One of the scientific goals of DXS and UDS is to 
locate Extremely Red Objects (EROs). Again, with 
a tiny fraction of the eventual survey data we are al- 
ready competitive with all previous studies. Using the 
same UKIDSS SV data as above, we cross-match with 
publicly available optical data from the INT Wide 
Field Camera Survejij, and select EROs as those with 
R-K>5 down to a limit of K=19. This produces a 
sample of 1660 EROs. The reliability of this sample 
is limited by the optical data at R=24, not by the 
IR data at K=19. Fig. 1131 shows the number counts 
of EROs in these data. Again, we have already dupli- 
cated most previous work in these science verification 
data. 



7 SURVEY PROGRESS 

The final version of this paper is being written in 
April 2007. We have made an 'Early Data Release 
(EDR)' (Dye et al. 2006), and First and Second Data 
Releases (DRl and DR2; Warren et al. 2007a,b). The 
EDR and DRl papers contain considerable detail on 
the areas covered, data quality achieved, QC cuts 
applied, example SQL queries for extracting science 
ready data, and so on, with the DR2 paper providing 
an incremental update. DRl contains approximately 
7% of the expected final dataset, with some surveys 
proceeding faster than others. Table[6]summarises the 
content of DRl. 



7.1 Survey quality achieved 

After filtering data through the QC process, UKIDSS 
is achieving its design goals, with the partial excep- 
tion of depth in stacked surveys, which is expected 
to improve in later releases. Except where noted, the 
statistics quoted below refer to measurements made 
from the DRl database. Distributions of these quan- 
tities are shown in Warren et al. (2007a). 

Image Quality. The median seeing in DRl 
is 0.82tn'i^", and the median stellar ellipticity is 



0.07 



+0.04 



(Icr ranges). These figures are expected to 
improve slightly in later releases ; nights used in the 
first half of DRl had imperfect alignment of the op- 
tics and fiexure correction. In QC filtering, seeing cuts 
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Figure 12. K-hand number counts from six hours of 
UKIDSS DXS/UDS SV observations of a single tile in 
the EL A IS Nl field. Results from the literature are shown 
for comparison. Note that the UKIDSS number counts in- 
clude both stars and galaxies. 
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Figure 13. Number counts of Extremely Red Objects 
(EROs), taken from the same DXS/UDS SV field as in 
Figure [TM but selected to have R — K > 5.0. These are 
compared with ERO counts from the literature, although 
note that Smith et al. (2002) use a slightly different se- 
lection criterion (R — K > 5.3). 



were applied to the various surveys, roughly as indi- 
cated in Tabled 

Photometric Accuracy. Calibration uses 
2MASS stars over the whole of each frame, and the 
internal dispersion indicates that locally we are reach- 
ing the target accuracy of 0.02 mag RMS., and that 
frame gradients are negligible. We also have compared 



the average colours of stars from field to field, which 
gives colour accuracies of 0.02 - 0.03 mag, indicat- 
ing that photometric errors are less than 0.02 mag in 
all filters. The global uniformity is currently as good 
as the 2MASS survey itself allows, which is believed 
to be 2% (Nikolaev et al. 2000). As well as cahbrat- 
ing from 2MASS, we are collecting observations of 
UKIRT faint standards, with which we will check uni- 
formity at faint magnitudes (Hodgkin et al. in prepa- 
ration). All the above applies to an average star in an 
average field. The analysis of DRl found some small 
but clear variation depending on stellar type and lo- 
cal extinction, indicating that some improvement is 
needed in colour transformations between the 2MASS 
and UKIDSS systems. These will be applied in DR2. 

Astrometric Accuracy. Astrometric accuracy 
was considered in some detail in the EDR (Dye et al. 
2006). As with photometry, the astrometric calibra- 
tion is tied to 2MASS stars, which in turn are tied 
to the Tycho system. The global uniformity there- 
fore assumes that of 2MASS, although this will in 
due course be checked independently. The precision 
of the astrometric solution changes with Galactic lati- 
tude, because of the number of 2MASS stars available 
for solution. The RMS accuracy varies from 50mas at 
low latitudes to lOOmas at high latitudes. 

Depth. The depth achieved in the shallow sur- 
veys is close to expectation. From LAS data in 
DRl, the median 5a point source depth in the stan- 
dard 40 second shots, compared to target depths is 
Y=20.16 (20.3), J=19.56 (19.5), H=18.81 (18.6), and 
K=18.19 (18.2). From the GCS, the Z depth achieved 
is Z=20.36 (20.4). For most surveys, this is a good 
guide to actual achieved catalogue depth, but in due 
course this will be analysed more fully in terms of 
completeness and relaibility. For the GPS in particu- 
lar, the effective depth is very sensitive to confusion, 
and so varies with position on the sky (see Warren 
et al. 2007). For the stacked surveys, experience so 
far has borne out what we found in the SV data - 
the depth achieved falls somewhat short of the naive 
y/N extrapolation. In the current pipeline output for 
the DXS, we are a few tenths of a magnitude short of 
target. At the eventual depth of the UDS, K=23, we 
expect of the order 100 galaxies per sq. arcmin, which 
should produce significant correlated confusion noise, 
as seen in other very deep IR surveys, such as the 
ISAAC FIRES survey of Labbe et al. (2003). How- 
ever, we are still some way short of this depth. The 
reason for the stacked survey depth discrepancy is 
currently unclear. It is possible that pipeline improve- 
ments and more careful analysis will remove much of 
this discrepancy. Further analysis will be provided in 
later data release papers. 

Database utility. The VDFS pipeline and 
archive system has delivered more or less what was 
required for UKIDSS, including reduced data quality 
and the ability to perform catalogue queries online. 
There are known processing limitations, listed on the 
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Survey 


Area deg 

all (any) 

filters 


Filters 


K-depth 
5(7 Vega 


Comf 


LAS 


190 (475) 


YJHK 


18.2 


6% 


GCS 


52 (401) 


ZYJHK 


18.2 


13% 


GPS 


77 (362) 


JHKH2 


18.1 


7% 


DXS 


3.1 (6.9) 


JK 


20.7 


11% 


UDS 


0.8 (0.8) 


JK 


21.6 


4% 



Table 6. Summary contents of the UKIDSS First Data 
Release (Warren et al. 2007). Depth achieved is shown 
only for the K-band, and the 'all' area is that with data 
achieving that depth and covered in all relevant bands. 
The completion estimate takes account of all filters, and 
areas with partial depth. 



WSA web page. These are relatively minor, with two 
main exceptions. The first is that profile magnitudes 
have not yet been implemented, which is potentially 
important for GPS analysis. The second is that deep 
stacking is not yet optimal. For this reason, the UDS 
working group have created their own image stack 
and deep catalogue outside the VDFS system; how- 
ever this has been supplied to the WSA and is cur- 
rently used as the publicly available product. Both of 
these issues will be resolved for later releases. 

Reliability and completeness. We have so 
far made only a preliminary study of completeness 
as a function of depth, or the number of spurious 
sources, for the various UKIDSS products. This is 
partly because they are databases rather than sta- 
tistically defined catalogues. A number of different 
catalogues can be extracted from the databases, de- 
fined with different cuts on quality parameters etc. 
No single catalogue is perfect - rather these should 
be extracted by users of the UKIDSS databases for 
specific scientific purposes. However, as the UKIDSS 
consortium, we will provide the best possible quality 
fiags, and in due course we will produce an analysis 
of representative catalogues as a guidance to users. 
Given that many scientific uses of UKIDSS will in- 
volve looking for rare or anomalous objects, reduc- 
ing contamination by spurious sources due to data 
artefacts, and wrong colours due to incorrect band- 
merging, are both crucial. The frequency and nature 
of artefacts is discussed in Dye et al. (2006) and War- 
ren et al. (2007), and source merging is discussed in 
detail in Hambly et al. (in preparation), but we have 
not yet estimated the density of spurious sources or 
incorrect mergers. Source variability can also mimic 
unusual colours, given that some passbands are mea- 
sured at separate epochs. This can only be assessed by 
attempting the relevant science projects. When this 
issue clarifies, later data release papers will include 
an estimate of the problem. 



7.2 Survey completion rate 

Up to DRl release, the rate at which the survey 
has accumulated has been slower than hoped. Raw 
data has been collected on average at 79% of the 
rate planned, and in addition 20% of data frames 
have failed our QC filtering, so that the net survey 
progress rate has been 63% of our goal. Since comple- 
tion of DRl observing, two major problems - bright 
moon ghosts, and an unexpectedly high data acquisi- 
tion system overhead - have been solved. Our current 
estimate is that the net survey progress rate is now 
~80% of goal. Within the first two years of UKIDSS, 
we expect to have completed approximately 2/3 of 
our two year plan. To achieve the full design param- 
eters of Table [3] will take approximately 1150 nights. 
The fraction of UKIRT time dedicated to WFCAM is 
under review, so the number of years required, and/or 
the adjusted survey parameters, are currently uncer- 
tain. 



7.3 First utilisation of UKIDSS 

In section |S] we have shown analyses undertaken by 
the consortium with science verification data. Since 
EDR and DRl, the first open community exploitation 
of UKIDSS has begun. Use of the UKIDSS databases 
in the WSA has been extensive. As of the end of 2006, 
over a billion rows of catalogue data have been down- 
loaded, and interesting results are starting to emerge. 
The first very high redshift quasar (z=5.86) has been 
found by Venemans et al. (2007). McLure et al. (2006) 
have reported the discovery of nine of the most lumi- 
nous known Lyman break galaxies, at 5 < 2: < 6 - 
they are surprisingly massive so soon (1.2 Gyr) after 
the Big Bang. A sample of 239 Distant red Galaxies 
(DRGs) at 2 < 2 < 3 has been selected by Foucaud 
et al. (2007) and their angular correlation function 
measured. Cirasuolo et al. (2007) have measured the 
evolution of the K-band galaxy luminosity function 
to z ~ 2. The coolest known brown dwarf (classi- 
fied T8.5) has been found by Warren et al. (2007c; 
spectrum shown in Warren et al. 2006). Lodieu et al. 
(2007) have found 129 new brown dwarfs, a significant 
fraction of the total known, including a dozen below 
20 Jupiter masses. Large new samples of quasars, dis- 
tant galaxy clusters, and extremely red objects have 
been constructed by Chin et al. (2007), Van Breukelen 
et al. (2006), and Simpson et al. (2006) respectively. 



7.4 Survey releases 

Data access poficy for UKIDSS is set by the UKIRT 
Board, and is set out on the JAG web page^J. 
UKIDSS is intended to produce multi-use data of 



'' http://www.jach.hawaii.edu/UKIRT/surveys 
/UKIDSSdatapohcies.html 



© 2007 RAS, MNRAS 379, 1599-1617 



The UKIRT Infrared Deep Sky Survey (UKIDSS) 23 



general benefit to astronomers worldwide, but with 
a temporary advantage for the communities that de- 
veloped the camera and surveys. Initially this meant 
UK astronomers, but now means any astronomer 
currently working in an ESO member state. The 
general principle is that the data are freely avail- 
able to any such Eistronomer from the point of re- 
lease, and available worldwide eighteen months later. 
(Note that individual members of the consortium 
have no privileged data access.) During the ESO- 
restricted phase, data access requires registration 
with the WSA. This is organised through a set of 
"community contacts" at astronomical institutions 
in ESO member states, who maintain their own 
databases of local users through the WSA system. 
Any reader who is not yet registered who believes 
they are eligible should contact their local commu- 
nity contact, or if necessary ask for a new community 
contact to be established. Fuller instructions and a 
list of current contacts is on the UKIDSS web page 
( |http ://www.ukidss.org/archive/archive.html). 

It is intended to make UKIDSS data available in 
a series of well defined staged releases. As well as in- 
volving incrementally more data, releases will usually 
involve reprocessing of data from previous releases, 
with updated correction of artefacts and so on. Each 
release will therefore be documented by a paper de- 
scribing the contents and limitations of that release. 
The first preliminary release, avalable from Feb 10th 
2006, had a relatively small amount of data (about 
1% of the expected total), and several known imper- 
fections in the data processing. This was therefore 
labelled an "Early Data Release (EDR)". It is de- 
scribed in more detail in Dye et al. (2006). The first 
full data release (DRl) took place on July 21st 2006. 
It contains 7% of the intended final data volume, and 
has substantially better data quality than the EDR 
in a variety of ways. It is described in Warren et 
al. (2007a). At the time of final revisions to this pa- 
per, the Second Data Release (DR2) has also recently 
taken place, with updated details in an online-only 
paper (Warren et al. 2007b). 

Further releases are likely to take place 
thereafter every six months or so. Raw WF- 
CAM data are available through the ESO archive 
system (http://www.eso.org) and through the 
CASU site (http://archive.ast.cam.ac.Uk/wfcam/j. 
All of the UKIDSS processed images and cata- 
logues are accessible and queryable through the 
web-based WFCAM Science Archive (WSA : 
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groups, plus a small number of other individuals. In 
addition to the bulk of the consortium, there are 
others we would like to thank. First and foremost, 
the UKIDSS enterprise would be impossible without 
the staff of UKIRT, the staff at UKATC who buiff 
WFCAM, and the staff of CASU and WFAU who 
built the data processing system. Next, we would like 
to note that much of UKIDSS, including its scien- 
tific ambition, design, data fiow concepts, and release 
plan, has been built upon the preceding ideas and 
high professional standards of the SDSS and 2MASS 
teams. Finally we would like to thank the anonymous 
referee, whose detailed and thoughtful comments im- 
proved this paper significantly. 
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